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We present a first-principles study of a coherent relationship between the optimized geometry 
and conductance of a three-aluminum-atom wire during its elongation process. Our simulation 
employs the most definite model including semi-infinite crystalline electrodes using the overbridging 
boundary-matching method [Phys. Rev. B 67, 1953f5 (2003)] extended to incorporate nonlocal 
pseudopotentials. The results that the conductance of the wire is ~ f Go and the conductance trace 
as a function of electrode spacing shows a convex downward curve before breaking are in agreement 
with experimental data. 

PACS numbers: 73.63.Nm, 68.65.La, 73.21. Hb 



Atomic-scale wires, which can be produced by scanning 
tunneling microscopy and break-junction techniques, are 
the objects of considerable attention pj, since the prop- 
erties of such minute systems with low dimensionality 
can be different from those of bulks; as the size of the 
structure becomes smaller than the electronic mean free 
path, the electron transport becomes ballistic and the 
conductance is dominated by the electronic structure. To 
date, there have been a great deal of geometrical for re- 
vealing the mechanical and electrical behavior of such 
minute wires. Among the most exciting discoveries was 
that a single strand of gold atoms suspended between 
electrodes manifests conductance quantization in steps 
of G =2e 2 /h @ as electrode spacing is increased, where 
e is the electron charge and h is Planck's constant. Fur- 
thermore, other studies showed that an atomic-scale wire 
made of aluminum atoms exhibits more complicated and 
even more interesting features; that is, the plateaus of 
conductance steps are not flat but have a positive slope 
before the wire breaks 0, Q . This somewhat counterin- 
tuitive result implies that the conductance increases even 
though the electrodes are pulled apart. 

On the theoretical side, there are many first-principles 
studies for examining electronic structure and transport 
properties of atomic-scale wires using tight-binding ap- 
roximation and /or the structureless jellium electrodes 
"BSEHHElHEi, some of which reported re- 
sults consistent with experimental data. For more strict 
theoretical predictions on details of electron conduction 
properties, it would be mandatory to completely elimi- 
nate artifacts and employ a realistic model where elec- 
trodes consist of truly semi-infinite crystals. Recently, 
the overbridging boundary-matching (OBM) formalism 
[r3.[l5j |. which is a novel first-principles treatment of elec- 
tron transport properties for nanoscale junctions sand- 
wiched between bona-fide semi-infinite crystalline elec- 
trodes, has been developed. 

In this Letter, we implement elaborate first-principles 
calculations with the incorporation of the OBM formal- 
ism to elucidate a close relationship between the relaxed 
geometrical structure and electronic conductance of the 
3-aluminum-atom wire during the elongation process us- 



ing a model in which the wire is suspended between two 
semi-infinite Al(001) bulks. To the best of our knowl- 
edge, this is the first theoretical investigation of this re- 
lationship where the optimized aluminum wire is con- 
nected to genuine semi-infinite aluminum crystalline elec- 
trodes which are gradually separated. We find that this 
3-aluminum-atom wire has a conductance of ~ f Go, and 
there is only one channel that gives a dominant contribu- 
tion to the electron conduction at the Fermi level. The 
conductance trace as a function of electrode spacing de- 
creases in the early stage of the elongation process and 
shows a minimum; in the case of shorter electrode spac- 
ing, a loop current (LC) is induced around the center 
atom of the bent wire so as to boost electron transmis- 
sion. With an increase in electrode spacing, the LC weak- 
ens and the conductance reduces. When electrode spac- 
ing is still further increased, the wire forms into a straight 
structure with a one-dimensional character, which con- 
tributes to the enhancement of electron transmission. As 
a result, the curve of the conductance trace is convex 
downward before the wire breaks, and subsequently the 
conduction channels between the two electrodes disap- 
pear. 

First, we determine the optimized atomic structure 
of the 3-aluminum-atom wire suspended between the 
A1(00I) electrodes during its elongation. Our first- 
principles simulations are based on the real-space finite- 
difference method ff6| . which enables us to determine 
both the self-consistent electronic ground state and the 
optimized geometry with a high degree of accuracy by 
virtue of the timesaving double-grid technique 1 1 71 and 
the direct minimization of the energy functional 18] . The 
norm-conserving pseudopotentials |19| of Troullier and 
Martins [2£j are used to describe the electron-ion inter- 
action. Exchange-correlation effects are treated by the 
local-density approximation |2l| in the density-functional 
theory. The nine-point finite-difference case, i.e., the case 
of Nf—4 in Eq. (A2) (here and hereafter, the equation 
number refers to that in Ref. [l^L is adopted for the 
derivative arising from the kinetic-energy operator in the 
Kohn-Sham equation. We set a cutoff energy of 25 Ry, 
which corresponds to a grid spacing of 0.63 a.u., and fur- 
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FIG. 1: Schematic representation of the model. The rectangle 
enclosed with broken lines represents the supercell used to 
evaluate the optimized atomic configuration and electronic 
structure. 

ther take a higher cutoff energy of 224 Ry in the vicinity 
of the nuclei with the augmentation of double-grid points 

m 

For optimization, a conventional supercell, which is in- 
dicated by the rectangle enclosed with broken lines in 
Fig. is employed. The size of the supercell is chosen 
as L x =L y =15.l2 a.u. and L z =22.68 a.u. + 2d. Here, the 
z axis is taken along the wire axis, the x and y axes are 
perpendicular to the wire, L x , L y and L z are the lengths 
of the supercell in the x, y and z directions, respectively, 
and d is the average value for the projections of inter- 
atomic distances between adjacent atoms of the wire onto 
the z component. The wire is connected to the square 
bases with the side length of at both its ends, mod- 
eled after the [001] aluminum strands, and all of these 
components intervene between the electrodes produced 
from three atomic layers of the Al(001) surface. The dis- 
tance between the electrode surface and the basis as well 
as that between the basis and the edge atom of the wire 
is chosen to be 0.5ao, where ao (=7.56 a.u.) is the lattice 
constant of the aluminum crystal. The position of the 
central atom of the wire is optimized on the (110) plane 
while the other atoms are fixed during the first-principles 
structural optimization. We find that the wire exhibits a 
bent structure below electrode spacing L es of 26.46 a.u. 
When elongated up to L es =27.09 a.u., the wire manifests 
geometrical transition from a bent wire to a straight one, 
and finally distorts at L es ~28.35 a.u. |2^| 

We next explore the electron transmission of the 
wire at zero-bias limit with the Landauer formula 
G=tr(TtT)Go, where T is a transmission matrix. The 
global wavefunctions for infinitely extended states con- 
tinuing from one electrode side to the other are evalu- 
ated by the OBM method . For the determination of 
the effective potential required in the OBM formalism, 
the electronic structure of the transition region which 
is composed of wire, bases, and a couple of atomic lay- 
ers of electrodes (see Fig. ^ as well as the electronic 
structure of the aluminum crystalline electrode is calcu- 
lated self-consistently employing a standard supercell ge- 
ometry with a three-dimensional periodic boundary. The 
wavefunctions in the crystalline electrode region are com- 



FIG. 2: Channel transmissions at the Fermi level for the 
three-aluminum-atom wires as a function of electrode spac- 
ing. 



puted by the generalized eigenvalue equation of Eq. (19), 
and in order to circumvent a numerical instability caused 
by evanescent waves, the wavefunction-matching formula 
Eq. (11) at the interfaces between the electrode and tran- 
sition regions is described in terms of the ratio matrices 
i?'s, which are defined by Eqs. (25) and (B9). Further- 
more, for calculating the conduction properties of the 
aluminum wires with s- and p-valence orbitals, we here 
extend the OBM scheme to the case including the nonlo- 
cal parts of the pseudopotentials of Troullier and Martins 
[2fj| : since the nonlocal parts are represented by block 
matrices in the Kohn-Sham Hamiltonian when the real- 
space finite-difference method is applied, the nine-point 
finite-difference formula for the kinetic-energy operator 
is adopted so that the off-diagonal elements of the finite- 
difference formula can entirely cover the block matrices 
of the nonlocal parts. Consequently, the Hamiltonian 
is described by a block-tridiagonal matrix composed of 
4iV X j / -dimensional block matrices, where iV™ is the num- 
ber of grid points on the x-y grid plane |24| . 

Figure |21 shows the channel transmissions at the Fermi 
level as a function of L es . The eigenchannels are inves- 
tigated by diagonalizing the Hermitian matrix T^T 0. 
There are certain features that are common to the undis- 
torted 3-aluminum-atom wire: throughout all the ranges 
of electrode spacing below L es =28.35 a.u., the first chan- 
nel, which has the characters of s and p z for a straight 
wire, is widely open, while the second and third ones, 
which have the characters of p x and p y and are degen- 
erate for the straight wire, have small transmission ~ 
1.0 x 10 -3 , and the resultant electronic conductance is 
~ 1 Gq. Moreover, one can see that the conductance 
trace as a function of electrode spacing exhibits a convex 
downward curve having a minimum at L es =26.46 a.u. 
Our result regarding the conductance trace accords with 
the available experimental data UB], and in particular, 
is in excellent agreement with the data recently reported 
by Mizobata et al. 0]. 

To gain an insight into this unique conductance trace, 
we plot in Fig. [3] the electron-density distributions and 
current densities of the first channel at the Fermi level for 
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several electrode spacings. The states incident from the 
left electrode are depicted. In the case of a short electrode 
spacing such as L es =24.57 a.u., the channel transmission 
is ~ 1, and the first channel current forms an LC rotating 
around the center atom. Upon elongating electrode spac- 
ing, one finds that the LC weakens and consequently, the 
transmission exhibits a minimum. The electron trans- 
port through the wire is closely associated with the be- 
havior of the LC: in fact, as the center atom is displaced 
upward from the most stable position depicted in Fig. 
the LC becomes enhanced to increase the transmis- 
sion. When the wire is stretched up to electrode spacing 
of 27.09 a.u., where it forms into a straight wire from a 
bent wire, the transmission recovers ~ 0.9. At this stage, 
the LC completely disappears and the current is axially 
symmetric around the wire axis. In the straight wire, 
the one-dimensional character strengthens, therefore the 
first channel transmission at the Fermi level is larger than 
that of the bent wire of L es > 25.83 a.u. After the wire's 
distortion, there are no conduction channels between the 
two electrodes. 

One of the previous theoretical studies on the conduc- 
tance of the 3-aluminum-atom wire employing structure- 
less jellium electrodes observed a similar current density 
@. However, this simulation does not have sufficient ac- 
curacy to discuss strictly quantitative properties of elec- 
tron transport because it was implemented using an infe- 
rior pseudopotential made only of local parts and a very 
restricted model where the bond lengths of the wire are 
fixed to be 5.4 a.u. For a definite interpretation of the 
relationship between geometry and conductance of the 
wire, a highly transferable pseudopotential having non- 
local parts and structural optimization are indispensable. 
Our result that the conductance is ~ 1 Go is in accor- 
dance with that of the ab initio tight-binding calcula- 
tion incorporating the square bases made of aluminum 
atoms 13] . Compared with the other studies using nonlo- 
cal pseudopotentials and structureless jellium electrodes 
without any atomic layers or bases, most of the chan- 
nel transmissions and conductance in the present study 



are smaller than those of the other studies: Lang |5j and 
Kobayashi et al. @ reported that the conductance of 
the wire is ~ 1.5 Go, and moreover, Kobayashi et al. 
reported the transmissions of the second and third chan- 
nels are ~ 0.3. In these calculations adopting jellium 
electrodes without any atomic layers, the channel trans- 
missions and electronic conductance sensitively depend 
on direct effects from the jellium electrodes, and there- 
fore the difference between our results and theirs may 
be mainly attributed to the absence of atomic layers or 
bases. 

In summary, we have studied a close relationship be- 
tween the optimized geometry and electronic conduc- 
tance of the 3-aluminum-atom wire suspended between 
genuine semi-infinite aluminum crystalline electrodes. 
Our findings are that (1) the 3-aluminum-atom wire 
sandwiched between two Al(001) electrodes exhibits an 
electronic conductance of ~ 1 Go, (2) there is one con- 
duction channel which widely opens at the Fermi level, 
and (3) the conductance trace as a function of electrode 
spacing behaves so as to draw a convex downward curve 
before the wire breaks. In the bent wire with a short elec- 
trode spacing of L es <26.46 a.u., the electron transport 
property of the wire depends on the strength of the LC. 
On the other hand, in the case of a long electrode spacing 
of L es >27.09 a.u., the one-dimensional character of the 
straight wire contributes towards the electron conduction 
to give rise to a conductance of ~ 1 Go- This result of 
the conductance trace is in accordance with experimen- 
tal data. Work in progress is the study on the electron 
conduction of the aluminum wire attached to the elec- 
trodes having other crystal faces such as the (111) and 
(110) planes. 
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